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A CASE OF ITERATION IN SEVERAL VARIABLES. 

By Albert A. Benhett. 

The Relation between Difference Equations and Iteration. 
1. If Fi(x) be a given function, the equation 

(1) /(* + 1) = FiVlx)] 

is a functional equation for fix), of the form of a difference equation, of 
the first order. The general difference equation of the first order is of the 
form 

f(x + 1) = H\fix), x], 

where H is a given function of two variables. The difference equation 

(1) has in any solution fix) obviously one arbitrary constant at least, 
since we may choose for the value of /(x ), say, y ; the values of /(x + n), 
where n is any positive (or negative) integer, will then be determinate. 
There remain, however, an infinite number of parameters in/(x), since the 
values of /(x + n) for n any real number between zero and one are still 
arbitrary. If we have any solution fix), we shall have 

fix + 2) -JMWCx)]} -F,[/(x)] 
and if we define jP„(x) for n, a positive integer, by the relation 

F n ix) m F^iFyix)] m F^F^x)), 

we shall then have for n a positive integer, 

(2) fix + n)= F n [fix)]. 

We may use (2) for the purpose of defining F n ix) for n other than a positive 
integer, if once fix + n) be defined for all real values of n, and the same 
feature of arbitrary definition f or < n < 1 appears in jF„(x). 

Conversely, starting with jP 1 (x), we may determine F„(x) for n a positive 
or negative integer and define it arbitrarily for every value of n between 
zero and one. The real solutions of 

fix + n) = F n [fix)], 

if Fi(x) be a real function defined for all real values of x, will contain 
but a single arbitrary constant. The problem of defining or determining 
F„(x) from Fi(x) is said to be a problem of iteration. Thus, for the case 
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of a single variable, the study of iteration is essentially the study of a 
special type of difference equation of the first order. 

We shall find that in the case of several variables, by a mere increase 
in the number of variables, we may make any system of difference equa- 
tions each of finite order, equivalent to a problem in iteration. Let us 
consider a set of p difference equations of orders ?i, q i} • • •, q p , respectively. 
Let one of these equations be of the form 

f(x + q) - H[f(x + q-l),f(x + q-2), • • ., f(x + 1), /(*) ; *]. 

By a change of notation, we may put 

x = f w (x), 

fix) =/(!)(*), 

/(* + l) -/(.,(*), 

f(x + q - 1) =/ ((t )(a0, 
so that we have the system 

/«,)(* + 1) = l+/o(z), 
/(!,(* +1) -/.(*), 
/<«>(* + 1) =/s(x), 

/(«-«(* + 1) -/«(*), 

/(«>(* + 1) = H[f iq) (x), /(«-«(*), • • •, /«(*), fi(x); /„(*)]. 

In this system, we note that x itself no longer appears explicitly. By 
repeating this process for each of the p equations, and by using through- 
out a single notation for x, viz., /(o>(x), we obtain a system of 1 + ?i 
+ ?2 ■+■ • • • + q P difference equations in the same number of unknown 
functions/, and each of the first order. Thus far, we have implied by our 
notation, that the number of variables occurring in any of the functions 
/ is one, and that these variables consist simply of x itself. If there be k 
variables, Xi, x-t, • • • , Xu, appearing in the unknown functions /, and in 
the functions H, we shall proceed as above and write 

Xi = /(o, o(xi, Xi, ■ ■ •, x k ), i = 1, 2, • • ■, k, 

and so obtain a system of k + q x + q 2 + • • • + q p = m difference equa- 
tions of the first order, in which the x's occur only within the functions /. 
After this normalization, there will be m equations in m unknown functions, 
each equation being of the first order, while the unknown functions will 
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contain k variables, xi, x i} • • • x*. We may next suppose that k — m. 
For if k < m, we may introduce as additional parameters into our un- 
known functions, m — k new variables, Xk+i, Xk+2, • • •, x m ; while if m < k, 
we introduce new equations of the first order in new unknown functions 
of the same k variables, and so obtain a system for which the number of 
functions and the number of equations are equal to the number of vari- 
ables. The reduction may be performed with respect to one variable at a 
time. Any system of difference equations of arbitrary orders in which 
the dependent variables for the greatest values of the arguments appear 
explicitly, together with certain new parameters, and functions, which 
may be taken in some simple form, may be reduced to a system of the 
following type: 



(3) 



fw( x i + 1, x 2 + 1, 
— F mi [f w (xi, 

/(2)(Zl + 1, Xi + 1, 

= FmiUm(xi> 



•, x n + 1) 
•, a; m ),/(2)(xi, 
•, x m + 1) 



•> Xm), * * •,/(») fa, '•', X m )], 



•) X m ), •'•, f(m)(Xi, '", X m )], 



• ',f(.m){Xl, •••, X m )], 

, x m - The equations 



f(m)(Xi + 1, X 2 + 1, • • •, X m + 1) 

— F(m)l[fm(Xi, • • •, x m ),fw(xi, • • •, X m ), 

the identity being in terms of the variables x r , xi, • • 

(3) we shall write symbolically as 

(4) f(x + 1) = F 1 [f(x)]. 

As in the case of one variable, we may regard the equation, now symbolic, 

/(* + n) - F n [f(x)], 

either as defining F„ in terms of f(x + n), or/(x + n) in terms of F„, the 
first being a problem of differences, and the second, of iteration. 

A Certain Class of Series and their Iteration. 

2. We shall investigate equations (3) in the case in which F^i[ui, u if 
• • • , u m ] is of the form 



(5) 



Fan[ui, u 2 , • • •, Um] = aui + (O 2 )), i = 1, 2, 



m, 



where by ((%*)) is meant a power series in the m variables (y,\, u 2 , • • •, u m ), 
for which no term appears of order less than k in the m variables taken 
together. In (5), the constant a is supposed to be the same for each 
value of i = 1, 2, • • •, m, and for the present we shall suppose that \a\ 
4= 0, * 1. 
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Let us suppose that we are given an arbitrary constant y, such that 
|y| # 0, and let us seek to determine the coefficients of m series, C, where 
C«) is of the form 

C ( i)[ui, u 2 , ••-, u m ] ss ym + ((u 2 )), i = 1, 2, • • •, m, 

such that the C"s satisfy the symbolic equation 

C[F 1 (u)] - Fi[C(m)]. 

The coefficients of C will be completely determined in terms of y, and of 
the coefficients of F t . For let us form a product of k u's distinct or differ- 
ent and call this product U. Then U will occur once and only once in 
each Fax with coefficients, respectively, say, d i} and once and only once 
in each C«) with coefficients, say, c,-. The coefficient of the term in U 
in C(i)[Fi(u)] will be of the form 

ydi + Cid h + (terms containing neither d{ nor c,-). 

The coefficient of the term in U in Fa)i[C(u)] will be of the form 

oxi + d,y h + (terms containing neither d t nor d). 

Furthermore c,- occurs for the first time in either C[Fi(w)] or F[C(u)] 
with the term in U, while all of the other terms in the coefficient of U, 
except that one containing d,- will contain exclusively coefficients of C 
and F, which appeared in these respective sets of series with terms in u's 
of degree less than h. Since we are supposing that the coefficients of F 
are known, and that, for the induction, the coefficients in C of terms in 
u's of degree less than h, are already determined, we have, to determine 
d, an equation of the form 

Ci(a h — a) = di(y h — y) + (known terms). 

Since \a\ 4= 0, H= 1, we see that c,-, also, is completely determined. The 
above argument does not apply to the initial term, and it is readily 
verified that y is indeed arbitrary provided that |?| #=0. 
Now F n (u) must be, itself, a special case of C(u), since 

F»[FM] = F 1 [F n (u)] = F n+1 (u). 

For n an integer, the initial terms of F n (x) are such that 

(6) F ii)n = a n Ui + ((w 2 )), i = 1, 2, • • -, m. 

We may choose for n other than an integer, say, 

F (i )n - a» +p ™ Ui + ((u 2 )), 
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where P(ri) is any periodic function of n of period unity and vanishing 
for n = 0. For the sake of simplicity, we shall choose for our definition 
of F(i) n , a series of the form (6) for all values of n, a choice which will be 
seen to be possible. In other words we shall choose for F n (u) the series 
C{u) for which y = a n . We shall need to use the fact that the form of 
F(i)n(u) is completely determined by this choice, and that the coefficient 
of any term in Fa)„(u) is a polynomial in a", while n enters in no other 
manner into the series. 

Application of Newton's Interpolation Formula. 

3. We shall now apply Newton's interpolation formula to the ex- 
pansion of F(i)„(u), and we shall prove that the series obtained from this 
formula, in a certain manner, about to be described, is not only, itself, 
a power series, but is identically equal to the form which would be ob- 
tained by the method of undetermined coefficients, as just considered. 

Newton's interpolation formula, for y = f(x), when the values 
Vi = f(xi) are known for i = 0, 1, 2, 3, • • •, is the following: 

f(x) = y + (x — x )[xqXi] + {x — x )(x — xJlxax&i] 
(7) 

+ (x — Xq)(x — Xi)(x — Xz)[XqXiX2X3] + • • •, 

where 

r i vi - yo r i y* ~ y* r 1 v* — y* 
fo* J = ^r^o ' lxiXi] = xT^x^ > M = x7=^ ' 

[xoXiXi] . [«*g = fr*J , [XlX2X3] = [«* J = M , 

r , [X\XzX^ — [XtiX\X^ 

[XtiXiXzXzi = x — x 

The only fact that we shall use concerning this formula is that it reduces 
identically to y„ if we put x = x n . To apply (7) to the question in hand, 
we shall put y = f(x) = Fa)„(u)/a n , x = a n , for every value of n, while 
the values of Xj, j = 0, 1, • ■ •, will be, respectively, a , a 1 , • • •, a>, 
The formula then becomes, 

^(»)n _ \ F(i) (a" — 1) / F(j)i _ F(i) \ 
a" ~ j 1 "*" (a - 1) \ a 1 / 

(8) 

(a" - l)(p» - a) ( Fun _ F un _ F un F im \ 1 

~t~ (a - l)(o 2 -1)U» a 2 o + 1 j 1 I 

or as it may be written, 
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(9) 



wAere V = F^o and 



(o n - 1) • ■ • (a" - a'- 1 ) 
+ (o _i) ... (o y_i) 



Fy+ •••], 



Vj{Ui, U 2) • • •, tt m ) = ^j Ky_i(Mi, w 2 , • • •, u m ). 

We shall first show that Vj(ui, u if • • •, w m ) is of the form ((w y+1 )). 
This is obviously true for j=0, and we shall assume it proved for j— 1 = 1, 
2, ♦ ♦ • , h. Let us take a term in U of degree h in the m's, and let the coef- 
ficient of U in Fa be 6. The terms obtained from U by replacing w,- 
by Fan, * = 1> 2, • • •, m in Vh, will, with one exception, be of degree higher 
than h. There will be a single term obtained from U of degree h of the 
form b(a h U). Hence in V h [F wu •••, F im) i]/a h+1 - V k (u u •••, u m ), 
the term in U will disappear, and since the same applies to every term of 
order h, we see that for j — 1 = h + 1, the above theorem as to the order 
of Vj is true. Hence every term of order h in the u's, which appears in 
the formula (9) for F«)„, will be obtained from V's of subscript less than h. 
The coefficient of any term of order h in (9) is therefore a polynomial in 
a" and is rational except for o" in the coefficients of F^i, F mu ♦•-, 
^(m)i. For n any positive integer (9) reduces to an identity and so must 
coincide with the series for F( i)n given by the method of undetermined 
coefficients. Since n enters in (9) only in the form a", the polynomial in 
a" which is in (9) the coefficient of a term in U will coincide for an in- 
finite number of distinct values of a", viz., o°, a 1 , a 2 , • • •, a 1 ', • • •, with the 
corresponding polynomial obtained by the method of undetermined coef- 
ficients. These polynomials must therefore coincide, and the expressions 
obtained by the two different methods must be identical. 

The Solution of the Difference Equation in Terms of the Iteration Series. 

4. Having obtained a form for F n , we may now secure a set of solutions 
/(,) of (3). For \a\ < 1, we shall have a definite set of series as the limit of 
a~ n F n , as n becomes positively infinite, while for \a\ > 1, the same formal 
series are obtained as the limit of a n F n f or n = + <» . The series in each 
case are 

GcoGfa, u>, • • ; «-) - [ V - ^~- V x + {a _ [^ _ 1} V 2 

(10) 

i - a - a 1 

F3+--J, 



(a - l)(a" - l)(o 3 - 1) 
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as is seen immediately. The G's satisfy the functional equations 

(11) aG(o(t«i, • • •, O = G (i) [F wl , F mu • • •, F (mn ], 
emce, for example, 

lim a(cr n F n ) = Urn cr n F n (F) 

and similarly for limit, n = — <» . Hence we have, 
1 + log„ G(o(w) = log„ (? (0 [iPJ 
or, replacing logo <?(o( w ) by a new notation, g^iu), 

(12) 1 + 0<o («) = ^oW 

We shall be obhged to restrict the set g^iu) to the extent that we shall 
suppose the system 

Xl = 0(1) (Wl, W 2 , • * •, Win), 

(13) x 2 = 0(2) (t*i, u 2 , •••, u m ), 

X m = g{m){Ui, u 2 , •••, u m ), 
to possess an inverse, which we shall denote by 

Ml = 0(1) -1 (^1, X 2 , •••, X m ), 

(14) U 2 = 0(2) -1 (*l> «2, • • *, Xm), 

U m = 0(m) _1 (^l> #2, • * ', X m ). 

Combining (12), (13), and (14), we have 

0(o -1 [l + Xi,l + X 2 , • ■ •, 1 + x m ] 

(15) _ 

= FmilSw 1 (xi, • • •, x m ), • • •, (m) \xi, ' • •, a: m )]. 

Hence, the set g^f 1 , i = 1, 2, • • •, m, constitute a solution of (3). 

Consideration of Some Limiting Cases. 

5. When o becomes equal to unity, the expression (9) assumes an in- 
determinate form, which may, however, be readily evaluated. A form 
for Fa) n when a = 1, may also be obtained by a direct method analogous 
to that used above. We shall have, for a = 1, 

F(,i)n = [ F(i) + J (F(i)i — F(i) ) 

( 16 ^ n(n - 1) 1 

H j-72 — (^"(02 — 2F(,-)i + ^(«)o) + - ■ • j 
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or as it may be written 

(17) F ( ,, B = [Fo + (^)F 1 +(^)F 2 +...+(^)F,+ ...], 



± Fit 



<»)<>• 



where V = i^oo and 

Vj = Fan ~ ( i ) Fu)j-i + ( 2 ) ^<»'^'- 2 "~ 

We shall not consider the case in which a = 0. If \a\ = 1, but & 4= 1, 
.;' = 1, 2, • • •, our original formulae hold. If, however, \a\ = 1, and also 
a> = 1, for some positive integer j, a new case occurs. Let h be the small- 
est positive integer (different from zero), for which a h = 1. Then the 
first h terms of (8) and (9) present no difficulty. The (h + l)st term 
contains however in the denominator the factor (a h — 1), while this factor 
will also occur in the numerator of this term, for n any positive integer. 
If we break off (8) after the hth term the formula will continue to be an 
identity for n — 0, 1, 2, • • •, h — 1. If F^h = ^<oo, in which case we 
always have a h = 1, then the first h terms of (8) will continue to represent 
F(On for every integer value of n, and will therefore, by the general ar- 
gument, represent a possible definition of F( i)n for all values of n. If 
F(Dh #= ^(oo, we cannot define F(i) n for all values of n as a series of the 
type we are considering.* 

Discussion of the Convergence of the Series. 

6. We shall now demonstrate that (8) is convergent when F^i is, 
i = 1, 2, '—, m. Let us suppose first that \a\ > 1. Let us suppose 
Fa)i dominated by a function F^i all of whose coefficients are real and 
positive. From the method of obtaining the coefficients of ^(02 from those 
of F(di, Fmi, F( 3 )i, • • •, F( m )i, we see that the coefficients of F (l) 2 will be 
real and positive and will dominate those of i ?, (»)2. Similarly for n a 
positive integer F(, )n will be dominated by F^n. Now, if F^n be de- 
termined by (8) as a convergent series for all values of n within a certain 
region, we see that since 5" is an increasing function of n for a > 1, 
F(i) n will be an increasing function of n for any given values of the inde- 
pendent variables, and in fact each term of (9) is an increasing function 
of n. We must show that for any set of functions, F W i, F( in , • • •, F( m )i, 
whose expansions in power series contain no term with negative coef- 
ficients, the V's, obtained as in (9), will have no negative coefficients. Let 
us suppose, for induction, that Vj-i has no negative terms, then we may 
write 

* Cf . the detailed discussion in the case of one variable, by the author, Annals of Math. 
vol. 17 (1915), p. 35 ff. 
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Vj(u x , u 2 , -•-, u m ) = ^{Vj-iiaui, ••-, au m ) + <p} — Vj-i(u 1} u 2 , • • •, u m ), 

where <p consists of partial derivatives of F,_i(awi, • • •, au m ) multiplied 
by products of [F( in (ui, u*, • • -, u m ) — au { ] for different or the same values 
of i, by simply expanding Fy_i(F (1) i, F Wh ••-, F (mn ) by Taylor's formula 
for the case of m variables. Now every term of <t> has a non-negative 
coefficient, and since a > 1, and F,-_i(wi, u 2 , • • •, «„) contains, for terms of 
lowest degree, those of degree j, we see that (l/a,i)Vj-i(aui, avh, • • •, au m ) 
— Vj-i(ui, U2, ' • ', u m ) contains no terms with negative coefficients. The 
coefficient in (9) of V } - is itself positive, since, as we have remarked, a n 
is an increasing function of n. It therefore suffices to establish the con- 
vergence of (9) for a set of dominant functions. 

Putting Fa)i = aui/(l — bs), where s = u x + w 2 + • • • + u m) we have 

(181 V ^'-1 •(! + o)(l + o + o 2 ) • • • (1 + o + • • • + a'-Vtt,- 

1 [1 - bs][l - 6(1 + a)s][l - &(1 + 5 + a*)s) ' ' 

■ ■ • [1 - 5(1 + 5 + ■ ■ • + ft- 1 )*] 

But the series ^ (t)n may be summed in this case in finite terms. We have, 
by inspection, a possible F<,)n in 

a n Ui 



(19) F (i)n = 



1 r l -5" ' 

1 — &t ^S 

1 — a 



which must coincide with the sum of the series given by (9). Since (18) 
converges when written as a power series, for 6(5 n — 1)|«| < (3 — 1), 
the set of series F wu dominated by ow,/(l — 6s), a > 1, will converge 
at least within this region. For the case \a\ < 1, we need merely invert 
the series to obtain o' = l/o such that \a'\ > 1, while a' takes the place 
of a in the inverted series. For a = 1, the expression (18) continues to 
hold while (19) becomes 

1 — nos 

which converges for |»|6|s| < 1. The convergence for \a\ = 1, but 
a h 4= 1, h = 1, 2, 3, • • •, depends upon other conditions as well as the 
convergence of F ( ,-)i as is seen in the case m = 1, discussed in the previous 
paper. 

Princeton University, 
October, 1915. 



